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1 INTRODUCTION 



The new generation of Cosmic Microwave Background 
(CMB) experiments, such as NASA's Wilkinson Microwave 
Anisotropy Probe ( WMAP) and a number of balloon-borne 
and ground-based instruments, are revolutionizing cosmol- 
ogy by providing precision estimates of a number of fun- 
damental cosmologic a l parameters (e.g. Bennett et all2003t 



I Teemark et all 120041: iMacTavish et al. 
l200fj). Theupcoming launch of ESA's 



(Bcrsancll i et alJll996tlThe Planck Collaborationll2005l) will 
provide an opportunity for even more precise parameter es- 
timates. 

In addition to the temperature anisotropy, the CMB 
is expected to be partially polarized due to Thom- 
son scattering of the anisotropic radiation field near 
the surface of last scattering. The first definitive de- 



tections of CMB i 


Dolarization have recently been made 


( Kovac et alJ 


2002 


iKosrut et al.l 12001 IReadhead et al. 


2004; Barkats et al. 


20051: iMontrov et alJl2005l: iPaee et al. 


2006). The angular power spectrum of polarized fluc- 



tuations can provide a wealth of addition al cosmologi- 
cal information (see e.g. iHu fc White! fl99^) . Perhaps the 
most tantalizing prospect is that primordial gravitational 
waves from the epoch of inflation will leave a distinct 
divergence-free (or 'B-mode') signature in the CMB po- 
larization that may be detectable by future experiments 
and cleanly separated from the curl-free (or '_B-mode') 
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ABSTRACT 

We report the results of pointed observations of the prototypical ultra-luminous in- 
frared galaxy (ULIRG) Arp 220 at 850 fj,m using the polarimeter on the SCUBA 
instrument on the James Clerk Maxwell Telescope. We find a Bayesian 99 per cent 
confidence upper limit on the polarized emission for Arp 220 of 1.54 per cent, averaged 
over the 15 arcsec beam-size. Arp 220 can serve as a proxy for other, more distant 
such galaxies. This upper limit constrains the magnetic field geometry in Arp 220 and 
also provides evidence that polarized ULIRGs will not be a major contaminant for 
next-generation cosmic microwave background polarization measurements. 

Key words: galaxies: individual: Arp 220 - galaxies: magnetic fields - cosmic mi- 
crowave background 

dominant polarization signal | Zaldarri aea fc SeliakI Il997l : 
iKamionkowski. Kosowskv. fc Stebbin3l^97F ! 

To realize this potential, careful control of systematic 
effects, including foreground emission, is essential. A num- 
ber of studies have characterized the potential of extragalac- 
tic s ources to contaminate CMB anisotropy measurements 
(e.g.lToffolatti et al.lll998l: ITeemark et al.ll200ol : iTucci et all 
120041: iToffolatti et alJl200aT In general, these authors con- 
clude that extragalactic contamination is most important at 
small angular scales (high multipole moments), that current 
estimates are not precisely constrained by available measure- 
ments, and that these sources are unlikely to pose a major 
difficulty to future CMB measurements. 

IScott fc White] lll999l) analy zed the data from early 
SCUBA surveys at 850 fim (e.g. |S m ail. Ivison, & Blain 1997; 
IHughes et all Il998t IBarger et al.l Il998t lEales et alJ Il999f> 

and concluded that the Planck mission may be confusion 
limited at frequencies of 350 GHz and higher, and that clus- 
tering of faint sources may be a measurable signal. 

iBorvs. Chapman, fc Scottl il999T> provided the first 

limit of the contribution of SC UBA point sour ces to 

the CMB. Several studies (e.g. lHaiman fc Knoxl l200d : 



20051 ISnereel etdl 

Planck mission 



iGonzalez-Nuevo et aljl2005l) have modeled how clustering 
of such sources may be manifest in the background over a 
wide range of wavelengths. 

The potential of polarized extragalactic sources to 
contaminate polarized CMB anisotropies has been much 
less well studied, particularly in the sub-mm region 
ide Zotti et al.ll999l:lTucci et all200lh . In this paper, we ad- 
dress this concern by investigating polarized sub-mm emis- 
sion from one particular object. 
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Arp 220 is a prototypical member of the class of ultra- 
luminous infrared galaxies (ULIRGs), with a far-IR lumi- 
nosity of approximately 1.6 x 1O 12 L l|Soifer et, al.1 11984 
iLisenfeld. Isaak. fc Hills) l2000l) . At roughly 70Mpc (H = 
70kmsec~ 1 Mpc~ 1 ) distance, it is the closest member of this 
class and one of the brightest galaxies in the local Universe. 
It is therefore well suited for studies that would be much 
more difficult for other, more distant ULIRGs. 

Although our initial motivation for investigating the 
sub-mm polarization of Arp 220 was as a means to esti- 
mate CMB contamination, such measurements of polariza- 
tion are interesting in their own right. Studying regions of 
polarized dust emission is important as a means of probing 
the ma gnetic field geometry responsi ble for aligning the dust 
grains jLazarian fc Finkbeineil200St) . Related issues include 
understanding the source of sub-mm emission, galactic su- 
perwinds, internal dynamics and dust physics. 



2 OBSERVATIONS 

We observed Arp 220 with the Submillimetre Common User 
Bolometer Array (SCUBA, Holland et al. 1999) at the James 
Clerk Maxwell Telescope (JCMT), Mauna Kea on 2000 Au- 
gust 23 and again on 2001 March 16. Conditions in the 
2000 run were favourable, with an 850 /J,m opacity of ~ 0.3 
(tcso ~ 0.07). The sky was much less opaque in the 2001 
run, where we enjoyed roughly a factor of 2 lower opac- 
ity. Obser vations were condu cted using the SCUBA po- 
larimeter jGreaves et al-llioOSt) which consists of a rotating 
quartz half-wave plate in front of a fixed wire grid analyzer, 
mounted externally on the SCUBA dewar. 

We chose to perform the polarization observations in 
'photometry' mode as opposed to the more common imag- 
ing mode, in order to achieve higher on-source efficiency. The 
angular size of Arp 220's infrared luminous core is smaller 
than the resolution of SCUBA at even its highest frequency 
channel, and thus we are not missing any flux by performing 
a photometry-mode observation. Because SCUBA employs 
a dichroic beam splitter, data are collected simultaneously 
for the 450 /im and 850 /im arrays. Although we are mainly 
interested in the central array bolometers (denoted 'C14' 
and 'H7', for the short and long wavelength arrays, respec- 
tively), data from the remaining array bolometers were also 
gathered and used as a monitor of the atmospheric emission. 

The observations required several levels of signal mod- 
ulation. The array was chopped in azimuth at 7.8125 Hz in 
order to remove common mode atmospheric signal in the 
source and reference positions. The chop throw was 90 arc 
seconds for the 2000 observations and 40 arc seconds for 
the 2001 observations. After four, 1 second integrations, the 
telescope was 'nodded' in azimuth to match the fast azimuth 
chop and another set of four, 1 second integrations were 
taken. The difference between the chopped signal in both 
nods removes instrumental eff ects which are depe ndent on 
the secondary mirror position (Z cmcov et al.ll2005l) . 

After these integrations, the polarimeter half-wave plate 
was moved to the next in a sequence of 16 rotational posi- 
tions. A full set of rotational positions was thus obtained, 
consisting of 128 seconds of observations, with an elapsed 
time of 280 seconds, including overheads. Arp 220 was ob- 
served at elevation angles ranging from approximately 30 to 



85 degrees above the horizon. From the two runs, the to- 
tal on-source observing time was approximately 3.4 hours, 
consisting of 95 full wave plate cycles. As we explain later, 
however, some data were not included in the analysis. 

In the 2000 observing run, several levels of systematic 
controls were performed. Two different sequences of wave 
plate positions were used: the first sequence was the stan- 
dard one, with sixteen 22.5° rotational steps, one after an- 
other; in the second method, the angular sequence was 0°- 
315° in steps of 45°, and then 22.5°-337.5° using the same 
step. Observations with the two sequences were interleaved. 
The different wave-plate strategies were performed in an at- 
tempt to detect or limit the contribution of atmospheric fluc- 
tuations to the polarization signal. Additional observations 
were performed in a similar manner using a different 850 /im 
bolometer ('G15') centred on Arp 220. After reducing these 
data and not finding any significant difference between the 
different approaches, we decided that for the second run the 
more straightforward approach of using the default wave- 
plate position order and H7 as the primary bolometer was 
better suited for observations of Arp 220. 

A polarization 'standard', DR 21 was also observed for 
15 total wave-plate cycles. This galactic star-forming re- 
gion ha s previously b e en ob served to be bright and po- 
larized. iGreaves etafl J1999T) found 2.34 ± 0.27 per cent 
polarization at 800 um w ith a position angle of 20° ± 3° . 
iMinchin fc Murray! jl994l) report 1.8±0.3 per cent polariza- 
tion at a position angle of 17° ± 4°, also at 800 fim. DR21's 
relatively high flux and likely lack of variability make this 
source convenient for cross-checking polarization measure- 
ments. 

Absolute flux calibration and determination of the in- 
strumental polarization (IP) were provided by observations 
of Uranus (13 wave-plate cycles) and Mars (12 wave-plate 
cycles). The Uranus observations were conducted with the 
same variety of systematic checks as the Arp 220 observa- 
tions. Pointing was checked roughly once per hour using the 
blazar 1611+343. Pointing corrections were typically a few 
arcseconds, except near transit - because the JCMT is an 
alt-az telescope, it has difficulty tracking sources such as 
Arp 220 that transit near the zenith. Arp 220 reaches an 
elevation of ~ 85° at the latitude of the JCMT. 



3 DATA ANALYSIS AND RESULTS 

The raw data were reduced with a combination of rou- 
tines from the SCUBA User Reduction Facility (SURF 
IJenness fc Lightfootll99'8T) and our own software. SURF was 
used to provide data free from atmospheric signal, while 
our own code was devoted to estimating the polarization 
strength. The roles of the codes are described below. 

3.1 Preliminary processing 

Using SURF, the data from the two nod positions were sub- 
tracted, flat-fielded and then corrected for extinction. We 
used the polynomial fits to the CSO 225 GHz opacity from 
the JCMT web page 1 . The opacity at 225 GHz was then 

1 See http://www.jach.hawaii.edu for more information. 
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co nverted to the 450 and 850 /im bands using the relations 
in I Archib ald et aT] l|2002 j). The residual sky background was 
subtracted l|,Ienness et al.lll998T) . as estimated from the me- 
dian of the twelve lowest noise bolometers on ring 3 of the 
array, and additionally a few anomalous 5-a or greater spikes 
were removed. 

Since we expect that the polarization will be weak, par- 
ticular care was taken to ensure the quality of the data that 
went into the analysis. In Fig. Q we plot the time-streams of 
the central bolometer in each array for both runs. 

3.2 Systematic error control 

Among the many sources of potential systematic error for 
polarization measurements, atmospheric transmission fluc- 
tuations, pointing errors or drifts, and spurious pick-up from 
the telescope environme nt are perhaps the most important 
jffildebrand et alJl200(tl . 

Contaminating pick-up was evaluated by processing the 
data from an off-source bolometer through the entire analy- 
sis pipeline. We have done this for the 'H9' bolometer. The 
polarized signal from this bolometer was consistent with 
zero, at the 1-a level, and had a similar noise level to the on- 
source bolometer analysis. We therefore conclude that po- 
larized contamination from the ground or from atmospheric 
emission (which we would expect to contaminate the on- and 
off-source bolometers at a similar level) does not contribute 
significantly to our on-source observations. 

We have examined the time series data for evidence of 
atmospheric transmission fluctuations. These would affect 
an observation of a bright point source differently than an 
off-source pixel, and may be expected to add noise to the 
measurement. If such fluctuations were significant and had 
somewhat shorter time-scale than the wave-plate rotation 
period, they would introduce a spurious signal that could 
be confused with intrinsic source polarization. One may also 
expect that a decrease in atmospheric transmission would be 
correlated with an increase in atmospheric emission (if due 
to a short time-scale fluctuation in opacity). In examining 
the time series data for the on- and off-source bolometers 
for our observations, we have found no evidence for such an 
effect. We are only sensitive to the spatial gradients in off- 
source emission, however, because only the differenced data 
after the fast chop are available from SCUBA. 

By examining the full time series of both the 450 and 
850 /im data, we have concluded that some observations were 
drifting somewhat off-source. A key signature of such a drift 
is a significant decrease in the amplitude of the observed sig- 
nal, and a greater change at 450 /im than at 850 /im due to 
the differences in beam-size. Such a decrease is evident near 
sample number 1000 in the right hand panel of Fig. 1. We 
have excluded from our final analysis the 7 full wave-plate 
cycles of data between pointing checks near where this point- 
ing drift occurred. We note that the telescope was pointed at 
very high elevation angle during this period, where tracking 
is most difficult. A comparison with and without these data 
included shows that our final polarization result changes by 
less than 1-a. Very small pointing drifts that may exist in 
the remaining data are not expected to affect the final result, 
since each wave-plate cycle has a gradient removed from the 
fit to Q and U. 

A period of increased noise near sample number 1000 in 



the left hand panel of Fig. 1 is evident. A comparison of data 
before and after sky subtraction shows that this is due to a 
period of increased sky noise. Excluding this section of data 
make less than a 0.3-cr change in the Arp 220 polarization 
fraction, and this section of data has been retained in our 
final result. In order to test the robustness of our results we 
have experimented with removing additional short sections 
of data. Aside from the 7 wave-plate cycles mentioned above, 
we have not found the results to be sensitive to the removal 
of other sections, beyond the expected increase in overall 
noise. 

Our observations of DR21 formally detect 850 /im po- 
larization at approximately the 4-<r leve l and are consistent 
fat th e 1-cr level) with those reported bv lMinchin fc Murray! 
(1994). It is not clear to us, however, that our off-source 
chop or sky subtraction is always in an emission-free region, 
so our claim for DR21 is one of broad consistency with pre- 
vious measurements, rather than of precise determination. 



3.3 Polarization analysis 

Using our own software, the average signal and estimated 
error for the four integrations at each wave-plate position 
for a full wave-plate cycle were fit (using singular value de- 
composition) to provide an estimate of the total flux and 
the degree and direction of polarization. A linear gradient 
in flux was removed from each waveplate cycle in an at- 
tempt to decrease the effect of residual pointing drifts. The 
results for all the wave-plate cycles were then combined after 
discarding several waveplate cycles due to pointing drift as 
described above. This fitting procedure essentially mirrors 
that of the observatory's available software, but allowed us 
to pipeline the analysis and rapidly compare a number of al- 
ternative reduction schemes before concluding that the pro- 
cessing described here was sufficient. Among the alternatives 
investigated (and eventually passed over) were fits to partial 
and multiple waveplate rotations, discarding waveplate ro- 
tation cycles based on a goodness-of-fit criterion, discarding 
waveplate rotations that produced outliers in Q and U, and 
an attempt at regressing fluctuations as measured with the 
450 /jm data. An overall instrumental polarization (IP) that 
is contributed mostly by the JCMT wind screen must also 
be subtracted before arriving at the final source polariza- 
tion estimate. We removed the IP estimated from previous 
observations of Mars and Uranus (assumed to be unpolar- 
ized) of 1.06 per cen t at a position angle of 161° given by 
iGreaves et ail J20031) . Through our observations of Mars, we 
derive our own estimate of the IP of .89 ± 0.23 per c ent at 
154° ± 7°, which is consistent with the lGreaves et all (120031) 
value. Our final results are relatively insensitive to the pre- 
cise value of the IP, because of the range of elevation angles 
during which the Arp 220 measurements were made. Switch- 
ing between the two estimates of IP above gives less than a 
1-a change in t he final results for S tokes q and u. We use 
the values from lGreaves et alJ i2003T) for our final results. 

The data reduction process was run on all the planetary 
observations to produce an estimate of the system calibra- 
tion. We adopt a value of 425 JyV -1 for the results pre- 
sented here. There is a systematic uncertainty in this value 
from atmospheric and pointing effects, which we estimate to 
be approximately ±5 per cent based on the distribution of 
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Arp 220 23/08/2000 data 



Arp 220 16/03/2001 data 




Figure 1. Timestreams for the extinction corrected, sky-subtracted on-source bolometer from the two runs. The top two panels show 
the 450 and 850 (im signal (in Volts). The third panel shows the ratio of the 2 signals, and should be constant under ideal conditions. The 
bottom plot shows the measured 225 GHz CSO r (black circles) and the polynomial fit (dotted line), together with the elevation angle 
of the observation (solid line). Vertical dashed lines indicate where pointing checks were performed. Left: The 2000 data are shown, with 
a solid black vertical line at sample 1664, denoting the transition from observations taken with H7 as the on-source bolometer to using 
bolometer G15. Right: The 2001 data highlights the problems associated with observing the target near the zenith. In particular, the data 
between samples ~ 750 and ~ 1250 demonstrate that the source probably drifted away from the bolometer centre around transit. These 
data were not included in the analysis. The opacity was much lower in the 2001 observations, and one can actually see the instrumental 
polarization in the 450 /im data. 



the calibration data. An overall change in the calibration, 
however, will not affect the polarization estimate. 

The results from the first and second observing runs are 
consistent, with most of the statistical weight coming from 
the second run. Our formal statistical result for Arp 220 at 
850 /J,m is 

I = 751 ± 38 mJy 
q = 0.00582 ± 0.00350 
u = 0.00426 ± 0.00448, 

where q and u are the normalized Stokes parameters, and I is 
the Stokes intensity. Our int ensity measurement is consistent 
wit h the values measure d bv lLisenfeld. Isaak. fc Hills! (|200ffl 
and lDunne et all i2000l) . A naive estimate for the polariza- 
tion amplitude and orientation can be produced from these 
measurements, using p 2 — q 2 + u 2 and tan(2#) = q/u. This 
yields a polarization fraction of p = 0.72 ± 0.39 per cent at 
a position angle of 18.1 ± 15.4 degrees. This value, however, 
considerably overestimates the true l evel of polar i zation due 
to the biasing effect of noise (see iRicd Il947t Scrkowski 
| l958l: IWardle fc Kronberlll974l : ISimmons fc Stewardll98.4 
Hilde brand et, al]l200flb To counter the effect of noise bias, 
we can produce an estimate of the posterior probability dis- 
tribution using a Bayesian framework and an assumed uni- 
form prior on the degree of polarization and the position 
angle. This procedure yields no firm polarization detection 



and gives our final result for the 99 per cent confidence upper 
limit on the polarization of Arp 220: 1.54 per cent. 

We note that pointing control is a likely systematic lim- 
itation to these or similar measurements. A systematic wan- 
der of 2.5 arcsec during the 12 minutes it takes to com- 
plete a wave-plate rotation is sufficient to induce a 1 per 
cent polarized signal. These concerns are more serious for 
the 450 fim measurements, which have a smaller beam-size. 
We therefore do not report a polarization result for the 
450 ti m Arp 220 measurem ents. The coming SCUBA-2 up- 
grade faolland et all 2003) will minimize the effects of point- 
ing drift and atmospheric fluctuations with a filled detector 
array and a much faster wave-plate rotation speed. 



4 DISCUSSION 

The lack of sub-mm polarization signal in Arp 220 is 
intriguing, but perhaps not very surprising. In our own 
Galaxy, polarized dust emission h as been observed in a 
variety of molecular clouds (e.g. iFiege fc Pudritj l2000l : 
Matt hews fc Wilson! Eool iHoude et alJ l2004h and also in 
pre-stellar cores ijWard- Thompso n et alJ I200 1 *) . Magnetic 
fields in these environments can result in aligned dust grains 
which emit radiation with an E-field preferentially aligned. 
The expected level of polarizati on is typically a few per cent 
jHildebrand fc Draeovanl Il995lb For extragalactic sources 
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Figure 2. Theoretical CMB angular power spectra with potential contributions from polarized dusty galaxies and with noise estimates 
from Planck. The solid lines show the expected CMB power spectrum from (top to bottom) temperature fluctuations, E-mode polar- 
ization, B-modc polarization with lensing, and B-mode without lensing. A standard cosmological model constrained by WMAP three 
year data has been assumed, along with a tensor-to-scalar ratio of 0.1. Also shown are the pixel polarized noise estimates from Planck at 
353 GHz (upper dashed line) and 143 GHz (lower dashed line). The dotted line shows an estimate of the contribution of polarized dusty 
galaxies at 350 GHz assuming all the sources are 1.5% polarized with random orientation angles, the sources are not clustered, and that 
sources above 100 mjy have been removed. The two dot dash curves show two estimates from de Zotti et al. (1999) of the contribution 
of dusty polarized galaxies at 143 GHz, assuming 2% galaxy polarization. 



iGreaves et alJ i2000t) have reported a detection of polarized 
850 /J,m emission from M82, via resolved imaging observa- 
tions with SCUBA. Because of the dust grain alignment 
mechanism, sub-mm polarized emission traces the magnetic 
field geometry averaged over the beam-size. Arp 220's small 
angular size means that the magnetic field would have to 
be aligned over a significant fraction of the source in or- 
der to produce detectable polarization, otherwise the ran- 
dom orientations from a variety of distinct regions sam- 
pled by the JCMT beam would tend to cancel out. Note 
that the average polarization o ver the entire M82 S CUBA 
mapping of M82 carried out by [Greaves et al] (120001) corre- 
sponds to only about 0.4 per cent. For Arp 220 high resolu- 
tion interferometric maps in CO (1-0) and dust continuum 
JScoville et aljll99ll) show that the emission is extremely 
concentrated in a dense core. Perhaps we can therefore con- 
clude that there does not exist a simple aligned magnetic 
field geometry in the core of Arp 220. The extreme dust en- 
vironment in Arp 220's core, however, makes it dangerous to 



draw conclusions based on dust properties in more benign 
envi ronments like that in our own Galaxy or even M82. 

IJones fc Klebel lll989T) detect weak near-IR polarization 
in Arp 220, which they interpret as due to a simple screen 
of aligned dust in front of a bright nucleus. The reported 
value is 0.54 ± 0.18 per cent p olarization at jf-band, with 
a posi tion angle of 58° ± 11°. ISiebenmorgen fc Efstathioul 
(2001) report a mid-infrared detection of polarization due to 
absorption of 3.1 ± 0.9 per cent at 14.3 /im with a position 
angle of 62° ±9°. Both of these studies also conclude that the 
dust grain alignment must be inefficient, otherwise the high 
optical depth would lead to a much higher degree of polariza- 
tion. It is unclear, however, if Arp 220's relatively low level 
of pol arization is typical, since ISiebenmorgen fc Efstathioul 
( 2001) report mid-IR polarization levels as high as 8 per cent 
in other ULIRGs. 

At different wavelengths one might expect the core 
polarization to reflect the geometry of the nuclear disk, 
the separation of the double nucleus and the orientation 
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of outflows. Th e gas disc of Arp 220 is at a pproximately 
PA = 45° (e.g. IScoville. Yun. fc Brvantlll997l) . while mod- 



els of the core radio, mm and sub- mm emission suggest 
2 compon ents separated by ~ 1 arcsec, with PA ~ 80- 



100° 


fe.g. |Baan & Haschickl 


119951: IScoville. Yun. & Brvant 


1997 


: Eckart & Downes! 2001 


iMundell. Ferruit & Pedlarrov 


2001 


). Obviously our upper limit sheds little light on 



whether the magnetic field structure is related to any of 
these features - this will await more stable measurements, 
as well as higher resolution studies. 

A naive estimate for the level of polarized fluctuations 
from extragalactic sources is the product of the level of 
fluctuations in flux a nd the average percentage polarization 
ide Zotti et alJll999h . A more detailed calculation can be 
done using the source counts for submillimetre galaxies. The 
shot noise from these sources results in a flat contribution 
to the CMB temperature or Stokes / Ces: 

/•Scut 

C{ = / S 2 (dN/dS)dS, (1) 
Jo 

where dN/dS are the differential source counts and Ct is the 
angular power spectrum. The contribution to the E- and B- 
mode power spectra are equal and given by 

Of = Cf =p 2 Cj/2, (2) 

for sources having fractional polarization p (see Tucci et 
al. 2004 and references therein). Assuming p — 1.5 per cent 
and following the counts estimate of Scott & White (1999) 
we show in Fig. 2 the level expected if individual sources 
can be removed at the 100 mjy level. We also show two es- 
timates from de Zotti et al. (1999) at 143 GHz, assuming 
each galaxy is 2% polarized. 

In producing CMB angular power spectrum estimates, 
one can remove or marginalize over pixels in the map with 
clearly detected point sources. The contamination of con- 
cern is then generally due to the brightest sources imme- 
diately below the detection threshold. The 100 mjy flux 
cut corresponds to the 4cr dete ction limit for the Planck 
353 G Hz channel all sky survey llThe Planck Collaboration! 
120051) . The Planck point source detection limits, however, 
can differ depending on the assumptions about the fore- 
grounds, the specific method adop ted, and whether data a t 
multiple frequencies are used fe.g.. lVielva et aL^OOllfeoOS^ . 
We note that raising the flux cut level to 200 mjy would re- 
sult in roughly a 10 per cent increase in the contribution to 
the angular power spectrum. 

Also shown is an overall view of the level of expected 
CMB temperature and polarization fluctuations using a 
standard cosmological model constrai ned by the WMAP re- 
sults, constructed using CMBFAST dSeliak fc Zaldarriaeal 
1996). The B-mode fluctuation curve (the potential contri- 
bution from primordial gravitational waves) was constructed 
assuming a tensor-to-scalar ratio of 0.1. 

The Planck 143 GHz and 353 GHz pixel noise contri- 
bution to the polarized angular power spectrum sensitiv- 
ity levels are also shown in Fig. 2. The standard error per 
l-mode is derived from the Pla nck instrument sensitivity 
l|The Planck Collaborat!onll2005l) using 



where of is the standard error per mode in the polarization 
power spectrum, f s y~y is the fraction of sky observed and Cf 
is the angula r power spectrum of the CMB polarization sig- 
nal (s ee, e.g.. lKnoxll995l : lKesden. Corrav. fc Kamionkowskil 
2002). Nf is the pixel noise contribution given by 

ATP f 47VS2 Z 2 °b 2 IA\ 

N l =/sky— — e > (4) 

where s is the effective sensitivity, r is the total integra- 
tion time, and a - ;, is the instrument beam gaussian width. 
The noise estimates here are simple ones and do not ac- 
count for the effects of non-uniform sky coverage, galactic 
cuts, detector 1/f noise, etc. Nevertheless, the Planck noise 
levels are more than an order of magnitude above the polar- 
ized dusty galaxy level and we therefore conclude that such 
galaxies are unlikely to contaminate CMB polarization mea- 
surements with Planck even in the higher CMB frequency 
bands. 

Contamination in the 350 GHz range may be an is- 
sue, however, for future CMB polarization measurements 
that attempt to reach the 0.01 tensor-to-scalar ratio level, 
particularly at the higher multipoles. Experiments designed 
to measure the lensing contribution to B-mode polarization 
may also need to carefully consider the potential impact of 
polarized extragalactic sources. 

Our current state of knowledge, however, is still incom- 
plete. The above analysis relies on a number of extrapola- 
tions and simplifying assumptions which may not be correct, 
including: the unknown source counts at the ~ 100 mjy level; 
precisely at what level individual sources can be removed; 
the real distribution in the degree of polarization of such ob- 
jects; and whether clustering, potentially involving partial 
alignments of polarization axes, might be significant. De- 
termining whether Arp 220's level of sub-mm polarization 
is actually typical will await measurements of more objects 
and greater control of systematic effects, which will be pos- 
sible with the new SCUBA-2 instrument. 
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